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Abstract. We report the discovery of the optical and near-infrared counterparts to GRB 001011. The GRB 001011 error 
box determined by Beppo-SAX was simultaneously imaged in the near-infrared by the 3.58-m New Technology Telescope 
and in the optical by the 1.54-m Danish Teles cope ~ 8 hr after the gamma-ray event. Here we implement the colour-colour 
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discrimination technique proposed by Rhoads (2001) and extend it using near-IR data as well. We present the results provided 
by an automatic colour-colour discrimination pipe-line developed to discern the different populations of objects present in the 
GRB 001011 error box. Our software revealed three candidates based on single-epoch images. Second-epoch observations 
carried out ~3.2 days after the burst revealed that the most likely candidate had faded, thus identifying it with the counterpart 
to the GRB. In deep R-band images obtained 7 months after the burst a faint (R=25.38±0.25) elongated object, presumably the 
host galaxy of GRB 001011, was detected at the position of the afterglow. The GRB 001011 afterglow is the first discovered 
with the assistance of colour-colour diagram techniques. We discuss the advantages of using this method and its application to 
error boxes determined by future missions. 



1. Introduction 

Gamma-ray Bursts (GRBs) are intense flashes of high energy 
photons that occur uniformly distributed on the sky. They were 
discovered in 1967 (see Bonnell & Klebesadel 1996 for a dis- 



cussion of the first GRB detections), but due to the lack of rapid 
and precise localisations, their emission at other wavelengths 



was not detected until 1997 (van Paradijs et al. 1997; Frail et 
al. 1997; Costa et al. 1997). Thus, for 30 years they were not 



localised at longer wavelengths (X-ray, UV, optical, IR, radio) 
and their distances could not be measured. The breakthrough 
that occurred in 1997 can be attributed to the advent of the 



Beppo-SAX X-ray satellite (Boella et al. [1997]), thanks to its 
fast (within a few hours) and precise (few arcmin error ra- 
dius) localisations. Since the deter mination of the redshift of 
GRB 970508 (Metzger et al. |1997| ) another 17 secure spectro- 
scopic redshifts have been determined to date (n ot considering 
the supernova SN1998bw; Galam a et al. 1998 ) , ranging from 
z = 0.43 (Vreeswijk et al. EOOlfc Hjorth et al. ( |Z0004 |2000b| ) 
to z = 4.50 (Andersen et al. |2000|). For an additional 8 GRBs, 



optical counterparts have been found, but with no conclusive 
published redshift determinations. 

It is now widely accepted that at least the long-duration 
GRBs originate at cosmological distances (with the excep- 
tion of GRB 980425) with isotropic equivalent energy releases 
ranging fr om 10 51 to 10 54 erg (see Van P aradijs , Kouveliotou 
& Wijers ( 2000| ) for a review). Frail et al. (2001 ) have recently 
claimed that considering corrections for possible beaming ef- 
fects the range of the inferred high-energy release is restricted 
to a narrower band around 10 51 erg. Current models invo ke the 
collapse o f a m assive star into a black hole (Woosley 1993 , 
Paczyhski 1998 ) or t he me rging of two compact objects (e.g. 
Lattimer & Schramm 1974 ). The intrinsic brightness of GRBs 
and their afterglows allows one to probe the nature of their dis- 
tant host galaxies and potentially the star-f ormation history of 
the early universe (Lamb & Reichart 2000). 

Since the first detection of an optical counterpart to a 
GRB, most searches have been based either on comparing 
a single epoch image to Digital Sky Survey (DSS) images 
in order to search for new objects or by comparing images 
taken through the same filter at different epochs in order to 
find fading, transient objects. However, many optical transients 
are fainter than the DSS limit at the time of the first optical 
follow-up observations, and in some cases the afterglow light 
curves show "plateaus" that could disgu ise their transient na- 
ture (GRB 00030 1C; Masetti et al. Eoool GRB 001007; Castro 



Ceron et al. 2001 ) o r even brightness incre ases ( GRB 970508; 
Castro-Tirado et al. |l998j; Pedersen et al. |l998j). Hence, it is 



clear that the identification process would benefit from alter- 
native identification techniques. One such alternative is using 
colour-colour selection techniques similar to those used for 
quasar selection for many years (e.g. Warren et al. 1991 and 
references therein). Before the Beppo-SAX era colour-colour 
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diagrams were applied to deep late-time CCD images with the 
purpose of detecting the presence of quiescent objects with 
anomalous colours within GRB error boxes. When the prob- 
ability of finding one of these objects in a small GRB error 
box was low (< 10~ 3 ) the object was considered potentially 
GRB related. Examples of such objects are extragalactic (e.g. 
quasars and AGNs; Vrba et al. 1995| , Luginbuhl et al. 1995) as 
well as Galactic (e.g. neutron stars; Sokolov et al. 1995, flare 
stars; Gorosabel & Castro-Tirado 
Hudec & Christian 



1995). 



1998 



white dwarfs; Motch, 



1990t novae; Zharykov, Kopylov & Sokolov 



The first reported attempt to exploit colour-colour diagrams 
for early afterglow identification of GRBs was carried out at 
the beginning of the Beppo-SAX era with optical images of 
GRB 9 701 1 1 taken 19 hours after the trigger (Gorosabel et al. 
1998^ . However, GRB 970111 did not show any detectable 



optical emission and the feasibility of colour-colour discrim- 
ination techniques was not demonstrated. Recently, Rhoads 
( 2001 ) has presented a detailed discussion of colour-colour 
discrimination techniques and the feasibility of using them to 
distinguish GRB power-law spectral energy distributions from 
curved (blac k-bod y) stellar spectra. In an a posteriori analysis, 
Simon et al ( 2001 ) show afterglow optical colours cluster in a 
precise position on the optical colour-colour diagrams. In this 
paper we present the result of implementing such principles by 
applying an automatic colour-colour discrimination software 
pipe-line to data taken for GRB 001011 just ~ 8 hours after 
the trigger. We demonstrate that the colour selection techniques 
can be successfully applied to identify GRB afterglows, includ- 
ing the use of near-IR data. 

In Sect. H we describe the observations of the GRB 00101 1 
error-circle. In Sect. || we detail a candidate-selection method 
based on the study of the optical/near-infrared colours of the 
objects in the field. We illustrate the technique by describing 
its application to GRB 00101 1, which enabled the discovery of 
its optical/near-infrared counterpart. Section |] shows the char- 
acteristics of the GRB 001011 afterglow and its likely host 
galaxy. Finally, in Sect. || and Sect. ^ we discuss and list the 
conclusions of our work. 



2. Observations 

GRB 001011 was detected on October 11.6631 UT 2000 by 
both the Wide Field Cameras (WFC) and the Gamma-Ray 
Burst Monitor (GRBM) on board the Beppo-SAX satellit e 
and localised with an accuracy of 5' (Gandolfi et al. 2000a). 



The position was later refined to 2' prior to our observations 
(a 2 ooo = 18' l 23 m 4.32 s , 6 2 ooo = -50°53'56"4; Gandolfi et al. 
ZOOObh. 



Optical (R-band) and near-infrared (J and Ks bands) ob- 
servations were carried out with the 1 .54-m Danish Telescope 
(1.54D) and the 3.58-m New Technology Telescope (NTT), 
both at ESO, La Silla, on October 11.9700-12.0361 UT, be- 
tween 7.37 and 8.95 hours after the burst. Optical and near- 
infrared comparison images were obtained with the same tele- 
scopes during the following week. Deep R-band images were 
obtained at the 1.54D in April 2001 and at the 8.2-m UT1 of 
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Table 1. Journal of observations of the GRB 00101 1 optical/near-IR counterpart. 



Telescope 


Date UT 


Seeing 


Filter 


Exp. Time 


Mag. 






(arcsec) 




(sec) 




NTT(+SOFI, LF) 


11.9700-11.9828/10/2000 


1.5 


Ks 


15x60 


17.56±0.05 


NTT(+SOFI, LF) 


11.9834-11.9959/10/2000 


0.9 


J 


15 x60 


w^eio.os 1 


NTT(+SOFI, LF) 


11.9964-12.0092/10/2000 


0.7 


Ks 


15x60 


17.58±0.04 f 


1.54D(+DFOSC) 


11.9934-12.0361/10/2000 


1.2 


R 


5 x600 


20.99±0.05 t 


1.54D(+DFOSC) 


14.0029-14.0521/10/2000 


1.0 


R 


5 x600 


23.74±0.23 


NTT(+SOFI, SF) 


17.9794-18.0138/10/2000 


1.0 


Ks 


30x60 


> 21.3* 


1.54D(+DFOSC) 


19.3710-19.4064/04/2001 


0.9 


R 


10 x200 


> 24.3* 


1.54D(+DFOSC) 


20.3647-20.4242/04/200 1 


0.8 


R 


14 x200 


> 24.6* 


1.54D(+DFOSC) 


22.3741-22.4149/04/2001 


1.1 


R 


6 x600 


> 24.5* 


VLT(+FORSl) 


20.1618-20.1905/05/2001 


0.9 


R 


8 x300 


25.38±0.25 


★ 3a upper limit. 



i Images used to construct the colour-colour diagram. 



ESO's Very Large Telescope (VLT) in May 2001, 6-7 months 
after the burst. Table [j] displays the observing log. 

The near-infrared NTT observations were made with the 
infrared spectrograph and imaging camera "Son OF Isaac" 
(SOFI) using both the Large Field mode (LF, October 11) and 
the Small Field mode (SF, October 17). The fields of view 
(FOVs) in SF and LF modes are 2.'4 x 2.'4 and 4.'9 x 4.'9, re- 
spectively. The 1 .54D observations were made with the Danish 
Faint Object Spectrograph and Camera (DFOSC) which pro- 
vides a FOV of 13'.7 x 13'7. Consequently, both the optical 
and near-infrared images taken a few hours after the GRB trig- 
ger covered the entire refined WFC error box. The last set 
of deep images were acquired with the FOcal Reducer and 
Spectrograph (FORS1) mounted at UT1 (Antu) of the VLT 

2. 1. Photometric Calibration 



The field was calibrated in the near-infrared using observations 
of the standard stars sj9178, sj9013, and sj9106 at different air- 
masses. These standards and thus the near-IR measurements 
reported in the present paper are based o n the JHKs photo- 
metric system introduced by Persson et al. ( 1998 ). The derived 
zero-point error was 0.03 mag in both the J and Ks bands. In 
the optical the field was calibrated using observations of th e 
Landolt field SA 107 at different airmasses (Landolt 11992b. 



to the USNO-A2.0 catalogue. Therefore, the calibration used 
in the present paper is offset by 0.17 mag relative to the pre- 
liminary mean zero-point based on the USNO-A2.0 cat alogue , 
which was used in the discovery report (Gorosabel et al. 2000). 



The zero-point error for the R-band was estimated to be 0.02 
mag. In Table || we give the RJKs magnitudes for 6 secondary 
standard stars in the GRB 001011 field (see Fig. [j]). The pho- 
tometry of these stars was performed using the DAOPHOT- 
II package (Stetson 1987, 1997). The RJKs zero-point errors 
were added in quadrature to the measurement error derived 
by the DAOPHOT-II package, giving the error in the mag- 
nitudes shown in Tables [j] and |2|. Stars #1, #5, and #6 cor- 
respond to the USNO-A2.0 stars U0375-35233918, U0375- 
35222053, and U0375-35221378, which are listed to have R- 
band magnitudes of 15.3, 16.1, and 16.3, respectively. Using 
our calibration we find R=15.64 ± 0.04, R= 16.24 ± 0.04 and 
R= 16.33 ± 0.04 (see Table ||). Based on these three stars we 
have calculated a mean R-band offset of 0.17 mag with respect 



3. Candidate selection using colour-colour 
diagrams 

The main principle behind colour selection of optical after- 
glows (OAs) is that they have power-law spectral energy distri- 
butions, which can be distinguished from curved thermal stellar 
spectra in colour-colour plots. This is a principle that has been 
used to select quasars for many years (Warren et al. 1991 and 
references therein). This fact can be used to carry out a first 
colour-based selection of OA candidates. These candidates can 
subsequently be checked by second epoch observations aimed 
at detecting variability. The advantage of having an identifica- 
tion scheme that is not based on the transient nature of OAs is 
that it allows fast follow-up spectroscopy or polarimetry, also 
for OAs that are fainter than the DSS limit. 

All afterglows observed so far as well as theore tical firebal l 
(e.g. Sari, Piran & Narayan[l99§ Meszaros |200i] P iran ^00l| ) 
or cannonball models (Dado, Dar & De Rujula 2001 ) show that 
GRB optical/infrared afterglow decays are reasonably well de- 
scribed by a spectral index independent of time (F v ~ v^ 13 , 
(3 not a function of time). Thus, the colours should remain ap- 
proximately constant with time, and the technique is applicable 
any time after the GRB. In the case that (3 is a function of time, 
then the afterglow decay is not achromatic. Thus, if the im- 
ages used to derive the colours are very separated in time and 
(3 varies strongly with time (extremely unlikely for GRB after- 
glows), the magnitudes can not be shifted easily to the same 
epoch assuming a power law decay F u ~ t~ a , with an unique 
achromatic value of a for all the bands. However, if the im- 
ages are contemporaneous (or at least quasi-simultaneous as 
in the case of GRB 00101 1) the error introduced by the decay 
epoch-sc alin g factor is negligible, even if (3 depends on time 
(see Sec. 3T for colour error estimates due to the epoch-scaling 
factor). So, the technique is not restricted to achromatic after- 
glow decays, if quasi-simultaneous data are used. Therefore, 
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Table 2. Secondary standards in the field of GRB 00101 1 





RA(J2000) 


Dec(J2000) 


Ks 


J 


R 


1 


18 


23:12.45 


-50:53:34.4 


13.91 ± 0.03 


14.33 ± 0.03 


15.64 ± 0.04 


2 


18 


23:11.34 


-50:53:57.7 


16.47 ± 0.06 


17.34 ± 0.05 


20.26 ± 0.06 


3 


18 


23:03.47 


-50:53:41.8 


15.46 ± 0.04 


16.17 ± 0.04 


17.59 ± 0.04 


4 


18 


23:01.65 


-50:54:17.0 


15.92 ± 0.04 


16.71 ±0.06 


18.38 ±0.05 


5 


18 


23:00.72 


-50:55:13.3 


14.90 ± 0.03 


15.28 ± 0.03 


16.24 ± 0.04 


6 


18 


23:00.04 


-50:54:36.4 


14.97 ± 0.03 


15.36 ± 0.03 


16.33 ± 0.04 



the method should be valid at any epoch after the gamma-ray 
burst at least until the emission of the host galaxy becomes 
dominant, usually weeks after the gamma-ray event. 

3. 1. Description and implementation of an 
identification pipe -line 

Due to the transient nature of OAs it is of foremost importance 
to have a fast identification method once a GRB localisation has 
been released. It is therefore very important to try to find an op- 
timal strategy for the identification of the OA. The most widely 
used strategy during the last three years has been to obtain an 
image, most often in the R-band, and then to look for new 
sources by comparing with DSS images (e.g., GRB 980425, 
GRB 980519, GRB 990123, GRB 990510, GRB 000301C, 
GRB 000926). In many cases this method may still be a fast 
and efficient way to localise the burst, but it is limited to the 
brightest OAs in fields of relatively little crowding. 

Colour-colour selection is an ideal complement to the com- 
parison with DSS images since it can be done with one-epoch 
observations and therefore allows fast follow-up spectroscopy 
or polarimetry. We have therefore designed a fast reduction 
pipe-line that constructs colour-colour diagrams of objects in 
GRB error boxes. Before the pipe-line can be used we need 
to obtain three nearly simultaneous images taken in the R, J 
and Ks filters. This requires the use of more than one telescope 
and a possibility for fast data transfer, but for major modern 
observatories such as e.g. ESO and with modern network facil- 
ities this is not a serious limitation. The three images are then 
WCS (World Coordinate System; e.g. Mink 1997) calibrated 
and fed to the pipe-line. The error introduced in the colours C 
by considering quasi-simultaneous images instead of simulta- 
neous ones is of the order of AC ~ where AT is the 
total time-span covered by the set of images, T is the mean de- 
lay of the observation with respect to the gamma-ray event and 
a is the instant power-law decay slope (F„ ~ t~ a ). In case the 
time-span covered in the observations AT is a non-negligible 
fraction of the delay T, then all the measurements should be 
shifted to the same epoch assuming a value of a. In the case of 
GRB 00101 1 the colour-colour diagram was constructed using 
R, J, and Ks-band images which overlap in time (indicated with 
a dagger in Table [[]). For a value of a — 1.33 (as was found in 
the optical, see Sect. Q) the errors in the colours are AC < 0.1 
mag. 

The pipe-line uses several selection criteria or filters to im- 
prove the identification of candidates. The first of these filters 



is to reject sources outside the GRB error box. Once the spa- 
tially coincident sources have been found, aperture photome- 
try is performed in the three images and a colour-colour dia- 
gram is constructed for all the objects inside the error box. The 
colours of the objects are dereddened using the Galac tic ex- 
tinction given by Schlegel, Finkbeiner & Davis (1998) in the 
direction of the GRB error box centre. In addition, the pipe- 
line creates a list of high-redshift OA candidates, consisting of 
those sources detected in the two reddest bands but absent in 
the blue one. 

For a quick analysis, the discrimination algorithm can be 
flagged to only use the relative colour differences among ob- 
jects in the GRB error box. We can derive the two-dimensional 
probability distribution of the uncalibrated colours in the field 
and construct the iso-density contour levels. Therefore, given 
an uncalibrated candidate we can determine the probability that 
this object belongs to the calculated probability distribution, 
i.e., whether its colours are typical of objects in the error box 
or not. This relative comparison largely removes the systematic 
effect of the Galactic reddening and also eliminates the influ- 
ence of a variable (and colour dependent) transparency of the 
atmosphere, encountered during adverse observing conditions. 
Therefore, the colour-colour discrimination technique does not 
strictly require absolute photometric calibration. The upper left 
panel of Fig. |^ shows an R— Ks versus J— Ks colour-colour di- 
agram for the objects inside the GRB 00101 1 error box (open 
diamonds). As can be seen, the colours of the GRB 001011 
counterpart (star) are inconsistent with almost all the other ob- 
jects within the error circle. 

A more careful analysis requires a rough photometric cali- 
bration. A good representation of the distribution of colours of 
field stars is given by the 2MASS and USNO-A2.0 catalogues 
(hereafter 2MASS+USNO), so the pipe-line can be flagged to 
use this catalogue to obtain a well constrained sequence of non- 
transient point sources near the GRB position on the sky. In 
the top right panel of Fig. || we show with dots and contour 
levels the position of 10097 sources from the 2MASS+USNO 
database. These 2MASS+USNO sources were selected to be 
inside a 10° x 10° window centred on the GRB 001011 er- 
ror box coordinates. The size of the window has been chosen 
to contain ~ 10 4 sources, thus giving a large colour sample 
for the construction of the contour levels in a reasonably short 
computing time. In order to eliminate the possible reddening 
gradient inside the window each 2MASS+USNO source has 
been dereddened by its individual reddening, rather than ap- 
plying a mean reddening correction to all sources. 
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Fig. 1. Upper figure: R-band image of the counterpart (between tick-marks) and the secondary standards in the field (see Table 
|J). The cross indicates the position of the refined WFC error circle centre. Central row of figures: from the left to the right side 
the R-band decay is displayed; 1 1.9934-12.0361 UT Oct 2000, 14.0029-14.0521 UT Oct 2000, 19.3710-22.4149 UT April 2001 
and 20.1618-20.1905 UT May 2001 (VLT detection of the host galaxy). Lower row of figures: the first three panels show the 
Ks-band brightness evolution; 1 1 .9700-1 1 .9828 UT, 1 1 .9964-12.0092 UT and 17.9794-18.0138 UT Oct 2000. The fourth figure 
displays the J-band detection on 1 1 .9834-1 1 .9959 UT Oct 2000. 
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The 2MASS+USNO catalogue trace can be used to deter- 
mine a preliminary colour calibration of the field. The colour 
calibration can simply be done by matching the median of the 
colours of the point-like sources in the field with the median of 
the 2MASS+USNO colours. By using the median rather than 
the average colour of the field objects, any bias introduced by 
objects with anomalous colours, which may be present in the 
field, is minimised. Only point-like sources are used to deter- 
mine the median of the colours, so the effect of the possible 
presence of extended objects (galaxy clusters or nebular struc- 
tures in general) is also eliminated. As seen in the upper right 
panel of Fig. |the median of the colour for the 2MASS+USNO 
catalogue is located at J— Ks=0.54 and R— Ks=1.33, respec- 
tively. The standard deviation in (J— Ks) and (R— Ks) is small, 
0.35 mag and 0.75 mag respectively. For a typical Beppo-SAX 
error box size (10 — 100 arcmin 2 ), at least 25 objects would 
be detected in all bands. Therefore, a colour zero-point for the 
field can be established with an accuracy of 0.07 mag in (J— Ks) 
and 0.15 mag in (R— Ks). This accuracy is fully sufficient for 
applying absolute colour selection criteria. 

Most of the 2MASS+USNO sources are located in a se- 
quence stretching from (J-Ks,R-Ks) = (0.1,0.5) to (1,3). As 



we will describe in Sect. 3.3 the colours of the majority of the 
2MASS+USNO sources can be reproduced by stellar spectral 
templates. 



3.2. The colours of OAs and the selection criterion 

In the upper right panel of Fig. [| we have over-plotted 
the colours of GRB 001011 (marked with a star) and the 
colours of OAs (open squares labelled from 1 to 7) with re- 
ported quasi-simultaneous R, J, and Ks-band magnitudes, i.e. 
GRB 971214 (#1: Diercks et al. ~ 



1998), GRB 980329 (#2: Fruchter 



1998|; Ramaprakash et al . 
T999] Palazzi et al. |l998|), 



GRB 980703 ( #3: Vr eeswijk et al. [1999]), GRB 9912 16 (#4 : 
Garnavich et al. |2000| ), GRB 00 030 1C (#5: Jensen et al. |200l| ), 
GRB 00092 6 (#6: Fynbo et al. [200 la| ) and GRB 010222 (#7: 
Masetti et al. 2001 ). Each GRB colour has been corrected for its 



Galac tic foreground extinction (Schlegel, Finkbeiner & Davis 
1998b. GRB 000131 was detected in the R and Ks-band but 



in the J-band only an upper limit was reported (Andersen et al. 
200C ). In order to visualise that the method is valid at least up to 
z — 4.5 we have over-plotted with an arrow on an open square 
(labelled with an 8) the colour constraints of GRB 000131. As 
can be seen in Fig. || the GRB 000131 colours are likely con- 
sistent with the colours of the other seven OAs (shaded area). 
Most of these OAs have R— Ks and J— Ks colours that place 
them well away from the locus of the stellar sequence. For a 
given R— Ks colour the OAs are 0.5-1.0 mag redder in J— Ks 
than stars. The mean dereddened colours for the seven OAs are 
(R - Ks) GRB = 3.3±0.9, (J-Ks) GRB = 1.6 ±0.3. The po- 
sitions of the known OAs in Fig. ^|do not show any correlation 
with their redshift. 

The spectral energy distributions of GRB afterglows are 
fairly well described by pure power-laws in the optical to near- 




infrared range, F v 



Using synthetic photometry we can 



3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 



Fig. 3. The figure shows the reddening of the R— Ks colour of 
a pure power-law as a function of the redshift due to Lyman-a 
blanketing along the line of sight to the observer. The blanket- 
ing effects for the two extreme cases of (3 — 0.6 (long-dashed 
line) and f3 = 1.5 (solid line) are plotted. Most GRB after- 
glows are intermediate between these cases. The reddening of 
the (R-Ks) colour due to the Lyman-a blanketing starts to be 
noticeable for redshifts z > 3.7 (see text). For z < 3.7 the 
blanketing effect is absent. 



therefore determine an analytic expression for the colour as 
function of /3, which is independent of redshift: 



R - Ks = 1.525 + $ x 1.284 
J — Ks = 1.035 + P x 0.540 

In the right upper panel of Fig. ^ we have over-plotted a 
solid line that corresponds to a pure afterglow Spectral Energy 
Distribution (SED) with ranging from 0.6 to 1.5, as expected 
for afterglows. All the traces and objects (including afterglows) 
plotted in Fig. ^| have been dereddened by their correspond- 
ing Galactic extinction (Schlegel, Finkbeiner & Davis 1998). 
As seen, the observed OAs are consistent with the predicted 
colours for power-law SEDs. 

When z > 4 the effect of Lyman-a blanketing becomes 
noticeable in the R— Ks colour, while the J— Ks colour is unaf- 
fected until z ~ 8. We have modelled the effect on the R Ks 
colour by imposing Lyman-a blanketing on a power-law SED 
with spectral slopes in the range 0.6 < (3 < 1.5 (M0ller & 
Jakobsen 1990; Madau 1995). The resulting spectrum was inte- 
grated for a range of redshifts, using the transmission function 
of an R-band filter as a weight function. The result is shown in 
Fig. ||. The effect of Lyman-a blanketing in the R-band has its 
onset at z — 3.7 and is essentially independent of the intrinsic 
spectral slope. At a redshift of z ~ 6.5, the attenuation reaches 

5 magnitudes, implying that GRB afterglows are unlikely to be 
detected in the R -band at redshifts much beyond 6 (see Lamb 

6 Reichart 2000 for a discussion on the detectability of very 
high redshift GRBs). 
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The effect of the Lyman-a blanketing on the colour-colour 
space can be visualised on the upper right panel of Fig. ^. As 
the Lyman-a forest shifts into the R-band sensitivity wave- 
length region the value of R— Ks increases moving the solid 
straight-line (pure power-law) upwards (dotted, dashed and 
dash-dotted line). For simplicity we have only over-plotted the 
loci of a pure power-law affected by blanketing when the af- 
terglow is located at z — 5.0 (dotted line), 5.5 (dashed line) 
and 6.0 (dash-dotted line). As can be seen in Fig. || the effect 
of the blanketing is to further separate the colours expected for 
power-law SEDs from the non-transient objects. 

In order to establish a list of OA candidates we use the fol- 
lowing colour selection criteria: R— Ks > 2.0, J— Ks > 1.2. 
These criteria define the rectangular shaded area of Fig. ^. This 
area contains the expected colours for power-law OAs as well 
as OA SEDs modified by Ly-a blanketing and at the same time 
has very little contamination from the stellar sequence. 

3.3. Location of other astrophysical objects in the 
diagram 

In order to improve the discrimination method, the different 
regions which the non-transient classes of objects occupy in 
the colour-colour diagram must be determined. To achieve this 
goal we either use synthetic photometry with spectral templates 
from the literature (stars, galaxies and quasars, following Wolf, 
Meisenheimer & Roser 2001, Wolf et al. 2001) or use colour 



pected inside the refined WFC error box is low (~ 10 ; based 



observations from the literature (asteroids). 

For stars we calculate R— Ks and J— Ks colours based on 



the spectral atlas of Pickles ( 1998 ), which contains spectra of 
131 stars with spectral types ranging from 05 to M8. It covers 
different luminosity classes, but concentrates on main sequence 
stars, and it also contains some spectra for high metalicities. 
The colours for these 131 stars are shown with filled circles in 
the lower left panel of Fig. ^| 

The expected colours for a population of galaxies were cal- 
culated using the galaxy spectral templates by Kinney et al. 
( 1996 ). This library consists of ten SEDs averaged from inte- 
grated spectra of local galaxies ranging in wavelength from 125 
nm to 1000 nm. The input spectra of quiescent galaxies were 
sorted by morphology beforehand into four templates and the 
starburst galaxies were sorted by colour into six groups yield- 
ing six more templates. Colours were calculated for all indi- 
vidual spectra redshifted from to 2.0 in steps of 0.01. The 
derived galaxy colours are displayed in the lower left panel of 
Fig. ^ with dashed lines. 

The quasar library is designed as a three-component model: 
We add a power-law continuum with an emission-l ine co ntour 
based on the template spectrum by Francis et al. ( 1991 ). The 
intensity of the emission-line contour was varied only globally, 
i.e. with no intensity dispersion among the lines. The slope of 



the power-law continuum F v 



,-0 



was varied in 15 steps of 



A/3 = 0.2 ranging from (3 = -0.8 to (3 = 2.0. The library 
was calculated for 151 redshifts ranging in steps of Az = 0.02 
from z = to z = 3. 

The number of quasars with redshift z > 3 having R<23.5 
(the limiting magnitude of the discovery DFOSC image) ex- 



on the contamination estimates given in Sec. 3.4 and in the 
fraction of z > 3 quasars present in the 10 th recompilation by 
Veron-Cetty & Veron 2001 ), therefore synthetic quasar colours 
have not been constructed for redshifts z > 3. The synthetic 
quasar colours are represented with a dotted line in the lower 
left panel of Fig. [| 

For satellites looking for GRBs in or near the ecliptic plane, 
such as HETE-II, it is also important to consider the contam- 
ination due to asteroids, especially as these objects do not ap- 
pear in the sky surveys that are usually used (e.g. DSS). Based 
on the asteroid colour recompilation by McCheyne, Eaton & 
Meadows ( |1985| ) we derived a mean colour of (R — Ks) = 
1.49 ± 0.28 and (J - Ks) = 0.45 ± 0.11 based on a sample 
of seven asteroids. An independent sample of albedos of nine 
asteroids reported by Gaffey et al. (1993) gives (R — Ks) = 
1.27 ± 0.23 (J - Ks) = 0.57 ± 0.27, consistent with the 
first sample. Thus, the location of asteroids in the (J— Ks) vs 
(R— Ks) colour-colour diagram is far from the region where we 
find GRB afterglows, ensuring that the colour-colour technique 
is effectively screening out asteroids. The loci of the above 
mentioned two asteroid samples is shown by the shaded back- 
ground region of the lower left panel of Fig. ^. 



3.4. Rejection of non-transient sources with similar 
colours as OAs 

Although there are many stars in the field, the stellar con- 
tamination is low. As can be seen in the upper right panel of 
Fig. U the GRB colour-colour space region (shaded rectangular 
area) is almost completely beyond the 3er contour level of the 
2MASS+USNO catalogue. 

The overlap with gal axy colours is more significant. 
According to Huang et al. ( 2001 ) we would expect ~4x 10 3 , 
10 4 and 2.4 x 10 4 galaxies per square degree with R<21, R<22 
and R<23, respectively. However, in the 21 < R < 23 magni- 
tude range only ~ 20% of the galaxies are compact under ob- 
serving con dition s of l'/O (Calar Alto Deep Imaging Survey, 
Wolf et al. 2001). Thus, we have designed a selection filter 
to permit identification and elimination of all extended objects 
from the colour-colour plot. To apply this filter, the Full Width 
Half Maximum (FWHM) distribution for all objects in the field 
is derived. Next, the mode of the distribution is calculated, and 
all the objects beyond ±3<r from the FWHM mode can be 
rejected. Based on the NTT deep field catalogue (Fontana et 
al. 200C) we estimate that ~ 45% of compact galaxies have 



colours consistent with GRBs (see the lower right panel of 
Fig. ||). Taking into account these estimates we would expect 
~ 1, 3 and 8 compact galaxies brighter than R= 21, R= 22 
and R= 23, respectively, in the 2 arcmin radius WFC error 
box. Thus, we can eliminate most of the galaxy contamination 
especially for bright OAs. This demonstrates that, besides the 
colour, the magnitude of the source is an important parame- 
ter that together with the colour shows the unusual nature of 
a source. However, we do not wish to exclude faint candidates 
since one of the important quests for future searches is to ex- 



8 



J. Gorosabel et al.: The Afterglow of GRB 00101 1 



tend the detection sensitivity to fainter magnitudes (Fynbo et 
al. 2001b). 

Quasars are not so easy to eliminate because they are point- 
like sources and their SEDs can be described by power-laws 
similar to those of afterglows. Fortunately (in this c ontext) , 
quasars are rare. Acc ording to Boyle, Jones & Sh anks ( 1991 ), 
Hartwick & Schade ( |1990| ) and Wolf et al. ( |1999| ) the number 
of quasars per square degree with R< 21, R< 22 and R< 23 
is ~ 80, 200 and 400, respectively. The number of quasars in- 
side the GRB 001011 error box would then be 0.3, 0.7 and 
1.4, depending on the depth of the image (R< 21, R< 22 and 
R< 23). As mentioned, GRB afterglows exhibit spectral in- 
dices ranging from 0.6 to 1.5. Quasars instead show values 
of (3 between and 1. Hence, afterglows tend to be redder than 
quasars. Therefore, the quasar contamination is mainly due to 
the reddest quasars. Only one third of the quasars show indices 
j3 > 0.6 (Francis et al. 1991). The pipe-line is able to exploit 
this property and has an option to reject objects with spectral 
indices smaller than 0.6. The spectral index of each object is 
calculated by least-square fitting a power-law SED to the R, J, 
and Ks-band magnitudes. 

Afterglows which are reddened by Lyman-a blanketing are 
shifted into a region of the colour-colour diagram which is 
occupied by L-type brown dwarfs. Although the density of 
brown dwarfs is l ow, of the order of on e per square degree 
(Kirkpatrick et al. ^999\ Fan et al. |2000| ) down to a Ks-band 
limiting magnitude of 20, this complicates the secure colour- 
colour identification of high redshift afterglows in large error 
boxes. 

Therefore, we consider that for a first selection of candi- 
dates the impact of contamination by non-transient sources is 
not severe, at least for GRB error boxes smaller than ~ 10 
arcmin 2 and possibly up to ~ 100 arcmin 2 . 

4. Results 

4. 1. Detection of the optical and near-infrared 
afterglow of GRB 001011 

In the upper left panel of Fig. || we show with open diamonds 
the R— Ks and J— Ks colours for all objects detected in all three 
bands in the GRB 00101 1 error box. Most of the sources lie in 
the sequence also seen for the 2MASS+USNO sources. When 
the selection method described above was applied, six candi- 
dates were selected, but three of them were rejected because 
of their non-stellar shape (using the selection filter described in 
Sect. ET4l 



Second epoch observations performed on Oct 14.0029- 
14.0521 UT at the 1.54D showed that one of the three candi- 
dates had almost disappeared (see Fig. |lj). The remaining two 
sources are compact, non-transient and have colours that are 
consistent with the OA colour-colour space locus. Both objects 
can be seen in the upper left panel of Fig. [| as the two closest 
open diamonds to the GRB 001011 colours (star) and consis- 
tent with the shaded area. The magnitudes of the two objects 
are R=21.04 and R=21.38, so the contamination is consistent 
with the expected number of compact galaxies and quasars with 
21<R<22 estimated in Sect. 3.4 (1+0.3 < number of compact 
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Fig. 4. The figure shows the afterglow history in the Ks, J, 
and R-bands. The two open squares show the Ks-band detec- 
tions and the diamond the J-band magnitude. The symbol size 
of these points has been enlarged in order to make easier the vi- 
sualisation of the error bars contained in the symbols. The solid 
line represents the o.r = 1.33 decay between our two R-band 
measurements (filled circles). The triangle shows the Ks-band 
upper limit imposed on 17.9794-18.0138 UT Oct 2000 by the 
NTT observations. The error bars along the horizontal axis rep- 
resent the exposure time of the observations. As can be seen 
the second Ks-band observation, the first R-band detection and 
the J-band measurement overlap in time, so they have been the 
basis of the colour-colour plot presented in Fig. || (see the ob- 
servations indicated with a dagger in Table [j]). The dashed line 
shows that the NTT upper limit is consistent with an achro- 
matic decay (q;r = oks)- 



galaxies + quasars < 3+0.7). The colours of both objects can 
be reproduced by synthetic traces of quasars or galaxies (see 
Fig. H lower left panel), so they are likely quasars or compact 
galaxies. 

An astrometric solution based on 50 USNO A2-0 reference 
stars yields a 20 oo = 18' l 23 m 04.56V 2 ooo = -50°54'15'/8 
(uncertainty l'.'O) for the OA, indicated in Fig. [l| These co- 
ordinates correspond to the Galactic coordinates I = 343.686°, 
b = —16.535° and a Galactic colour excess E(B— V)=0.10 



(Schlegel, Finkbeiner & Davis |1998| ). 

The position of the counterpart is fully consistent with the 
position of the refined 2' WFC radius error circle (Gandolfi 
et al. 2000b) and is located only 19V 5 from its centre. The 
counterpart turned out to be the object with the reddest J— Ks 
colour among the 90 point-like sources located within the GRB 
error box and detected in the three bands. The dereddened 
colours of GRB 001011 are (J-Ks) GRB 001011 = 1.62 ± 0.06, 
(R— Ks)grb 001011 = 3.17 ± 0.07, fully consistent with the 
dereddened colours of GRBs detected to date in the R, J, and 
Ks-bands and close to the trace of OA power-law SEDs (see 
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solid straight line of Fig. ^| right upper panel). The counter- 
part was point-like in the optical and in the near-IR, showing 
the same point-spread function as neighbouring stars. In or- 
der to improve the photometry given by the automatic pipe- 
line more accurate photometry was carried out by means of the 
DAOPHOT-II photometry package. The final photometry of the 
counterpart is shown in Table [I] and displayed in Fig. Q The de- 
cay indexes in the R and Ks bands are qr = 1.33 ± 0.11 (solid 
line of Fig. and q;ks = 0.24 ± 0.73, respectively. The large 
error in the determination of «k s is due to the short time-span 
between the two Ks-band detections. 

In order to measure a possible deviation from an achromatic 
decay, we have calculated the expected magnitude difference 
between the second and the first epoch Ks-band detections if 
qr = qks is assumed. The prediction (AKs = 0.12 mag) is 
only 1.5 g different from the measured magnitude difference 
(AKs= 0.02 ± 0.06), so an achromatic decay is an acceptable 
approximation for our two Ks-band detections. A power-law 
fit to the two Ks-band detections, fixing an achromatic decay 
index («Ks = «r = 1-33 ± 0.11, see dashed line of Fig. Q), 
yields Ks= 21.85 ± 0.35 on 17.9794-18.0138 UT Oct 2000. 
Thus, the NTT upper limit (Ks > 21.3, represented with a tri- 
angle in Fig. Q) is consistent with an achromatic fading. We 
conclude that our measurements are consistent with an achro- 
matic decay, although we can not exclude a more complex 
Ks-band lightcurve, with plateau phases (e.g. GRB 000301C, 
Rhoads & Fruchter pOOl^ GRB 010222, Masetti et al. |200l[ ) or 
re-brightenings (e.g. GRB 971214; Gorosabel et al. |1998b| ). 

Once the RJKs-band quasi-simultaneous magnitudes (see 
magnitudes indicated with a dagger in Table |IJ) have 
been dereddened by Galactic extinction in the direction of 
GRB 001011 (E(B-V)=0.10; Schlegel, Finkbeiner & Davis 



1998) and shifted in time (assuming an achromatic fading with 



a = 1.33) to the mean epoch of the second Ks-band detec- 
tion (12.0028 UT), a least-squares fit provides a spectral in- 
dex of j3 = 1.25 ± 0.05. The value of (3 has been derived 
without considering intrinsic extinction, which is unknown for 
GRB 001011. Therefore, (3 = 1.25 has to be considered as an 
upper limit to the actual OA spectral index. 

Given that we have measured a and constrained the value 
of (3, these values can be compared to the predictions given 
by several afterglow models. However, given the poor cover- 
age of the light curve and the absence of an estimate of the 
intrinsic extinction we stress that only weak conclusions on the 
afterglow physical properties can be deduced from our data. If 
we assume an unextincted afterglow then the upper limit on 
f3 would be close to the actual spectral index. According to 
the adiabatic expansion of a spherical afterglow in the slow- 
cooling regime (Sari, Piran & Narayan 1998| ), a=3/3/2 (for 
v < v c ) and a=(3/3-l)/2 (for v > v c ) is expected. Thus, the 
prediction of a given for v < v c is a — 1.88 ± 0.08, which 
is 4(j away from the measured value. However, the predicted 
value given for v > v c is a — 1.38 ± 0.08 fully consistent with 
the measured value of a = 1.33 ± 0.11. Besides, the derived 
value of the electron power-law index, p = 2/3 — 2.5 ± 0.1, is 
in agreement with those obtained for other afterglows. For the 
more general and common case of an extincted afterglow the 
value of (3 would be decreased and it could be compatible with 




Fig. 5. The figure shows a contour plot of the co-added R- 
band images taken in May 2001 with the VLT. The total expo- 
sure time is 2400s with a mean seeing of Of! 9. The cross shows 
the position of the optical counterpart, fully consistent with the 
object. The object is clearly elongated in the North-West direc- 
tion. The image has been smoothed with a 3 x 3 pixel boxcar 
filter. North is to the top and East to the left. 



the predictions given by the spherical adiabatic expansion for 
the case of v < v c . The measured value of (3 = 1.25 ± 0.05 is 
also explainable i n the context of the cannonball model (Dado, 
Dar & De Rujula |200l|). 



4.2. The host galaxy 

When all the R-band images taken in April 2001 at the 1.54D 
were co-added, the presence of a very faint extended object 
~ 3(7 above the background was noticeable. Further deep opti- 
cal observations carried out one month later with the VLT con- 
firmed the presence of an object (detected at 8(7 significance) 
with R= 25.38±0.25 measured in a 3'.'0 diameter circular aper- 
ture. The coordinates of the object are 012000 = 18 ,I 23 m 04.58 s , 
S 20Q0 = -50°54'16'.'0 (uncertainty 0"5), which is 0'.'28 from 
the OA position. Thus, the OA seems to be centred on the ex- 
tended object. This strongly suggests that the object is the host 
galaxy of GRB 001011. A contour plot of the object and the 
position of the OA can be seen in Fig. || The host galaxy is 
elongated in the North West-South East direction (PA = -45°). 



5. Discussion 

The colour-colour discrimination pipe-line presented in this pa- 
per makes use of the R, J, and Ks-band magnitudes. However, 
it can be applied using other bands as well. The best choice 
of the bands would be those where most of the curved stellar 
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spectra are as distant as possible from the afterglow power-law 
SEDs. A detailed study of the optimal configuration of near- 
IR + optical filters is however beyond the scope of the present 
paper. The method can also be generalised to a higher number 
of filters which would allow one to construct an N-dimensional 
colour-colour space. 

Obviously, the colour-colour selection is only useful for 
GRBs that have optical and/or near-IR afterglows. However, 
with 8m class telescopes the technique can easily be used for 
bursts as faint as R = 24. Hopefully this technique will there- 
fore facilitate the detection of optical and/or near-IR afterglows 
from a larger fraction of well localised GRBs than the ks 30% 



during the latest 3-4 years (Fynbo et al. 2001b) 



Another possible bias could affect the selection of high red- 
shift afterglows. As long as the Lyman-a forest does not enter 
into the bluest observed band the method does not introduce 
any redshift selection bias. At a redshift z > 3.7 the Lyman-a 
forest enters the R-band, and the spectrum is no longer well de- 
scribed by a single power-law SED. However, the effect would 
only make the OAs redder in R— Ks and thereby move the OAs 
further away from the non-transient sources. 

Therefore, as long as the object is detected in the R band 
~ 5cr above the sky background, the blanketing does not repre- 
sent any inconvenience to locate and distinguish objects in the 
colour-colour diagram. Conversely, it would help to discrim- 
inate a candidate from the non-transient sources in the GRB 
error box, at least up to a given upper limit in the redshift. For 
higher values of the redshift where the blanketing effect makes 
that the R-band detection is below ha, then the large photo- 
metric errors would be a difficulty to distinguish the OA from 
the rest of the objects. In this regime the blanketing starts to 
represent an inconvenience, rather than an advantage. For even 
larger redshift values the blanketing effect can be so severe that 
the afterglow is likely to remain undetected in the R-band. 

In the case of GRB 001011 the R-band detection is about 
2.5 magnitudes above the detection limit of the image, hence 
we could have tolerated a blanketing effect up to z < 5.8 (see 
Fig. [3]). We estimate that for the set of three images (1.54D 
R-band, NTT J and Ks-band) used in the present study, the 
Lyman-a blanketing helps the discrimination of the afterglow 
in the redshift range 3.7 < z < 5.3. For a redshift range 
5.3 < z < 5.8 the GRB 001011 R-band detection would be 
below 5(7 having R-band photometric errors above ~ 20%, so 
the efficiency of the pipe-line goes down. For redshifts z > 5.8 
the blanketing would make the GRB 001011 afterglow unde- 
tectable for our R-band images. For other limiting magnitudes 
and combinations of filters, the corresponding redshift ranges 
will be different. 

Obviously, the former discussion is only dealing with the 
impact of the blanketing on the colours of GRB 00101 1, and it 
has not considered the monochromatic cosmological dimming 
factor given by Dl(z)~ 2 X (1 + ^ — 1 , where Dl(z) is the 
luminosity distance. This effect would enlarge the photomet- 
ric errors of the high redshift afterglows an d hen ce complicate 
their discrimination (see Lamb & Reichart 2000 for a detailed 
discussion on the effect of the monochromatic cosmological 
dimming factor). 



The application of the method to very reddened afterglows 
which are bright in the near-IR but with no detectable emission 
in the optical is possible as long as a constraining lower limit on 
the R— Ks colour can be derived. This could happen under con- 
ditions of high extinction or with extremely high redshift bursts 
having 5.8 < z < 8. The pipe-line considers the potential case 
of highly extincted and/or redshifted afterglows creating a list 
of sources only detected in the near-IR (both J and Ks bands), 
and hence absent in the optical image. 

The efficiency of the pipe-line in the 5.8 < z < 8 redshift 
range depends on the constraint imposed by the R-band image. 
As can be seen in Fig. ^| for an R-band limiting magnitude of 
R > Ks+2 a combined detection in both the J and Ks bands 
would be able to distinguish the shaded area from the rest of 
the colour-colour space. The discrimination of undetected op- 
tical sources in the R < Ks+2 region is uncertain. For z > 8 
the Lyman-a blanketing also affects the J-band and a single 
analysis becomes much more complicated. 

Finally, there is a bias against OAs situated in bright host 
galaxies, both because the flux of the underlying galaxy may 
move the integrated flux of OA + galaxy away from the shaded 
region in Fig. || and because the galaxy will make the OA look 
extended. 



6. Conclusion 

We present a colour based selection pipe-line of OA candidates 
that only requires three quasi-contemporaneous images. The 
colour selection software is designed to be used in parallel with, 
or as input to, the normal technique, searching for transient 
sources through the comparison of observations from (at least) 
two epochs or with the DSS. Furthermore, it is a method that 
allows a fast identification of candidates for follow-up spec- 
troscopy also for OAs that are fainter than the DSS limit at 
the time of first optical/near-infrared observations. Thus, the 
technique is not meant to replace the normal procedure, but 
to complement it. The combination of colour-colour and vari- 
ability information could be a very powerful mean of doing 
automated OA discovery, with great potential for forthcoming 
missions like Swift. 

The technique has several advantages; it can be applied any 
time after the gamma-ray event, using a single set of images 
in different filters, and it does not strictly require absolute pho- 
tometric calibration. Another additional advantage is that the 
method is independent of the redshift, at least for redshifts with 
a negligible Lyman-a blanketing along the line of sight. 

In the case that the discrimination method is based on the 
R, J, and Ks magnitudes, it becomes independent of the red- 
shift for z < 3.7. For redshifts 3.7 < z < 5.3 the possible 
blanketing effect along the line of sight helps the colour dis- 
crimination. In the range 5.3 < z < 5.8 the blanketing in- 
troduces considerable errors in the R-band photometry and it 
starts to be an inconvenient. For 5.8 < z < 8 the efficiency 
is still valid as long as a constraining lower limit on the R— Ks 
colour can be derived (R > Ks+2). For extremely high redshift 
afterglows (z > 8) the colour-colour discrimination technique 
is uncertain. 
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This kind of analysis is very suitable for small GRB error 
boxes such as the ones reported by the NFI of Beppo-SAX or 
for the ones that HETE-II is expected eventually to determine, 
in which case the field contamination of quasars and compact 
galaxies is small. 

Using this method we discovered the GRB 001011 after- 
glow. The GRB 001011 optical afterglow evolution is consis- 
tent with a decay index of an = 1.33 ± 0.11. With no correc- 
tions for the intrinsic absorption, we derived a spectral index 
of (3 = 1.25 ± 0.05. Therefore, this value of (3 has to be con- 
sidered as an upper limit to the unextincted afterglow spectral 
index. If we assume negligible intrinsic absorption, the values 
of a and (3 are consistent with a spherical afterglow model with 
an electron energy index p = 2.5. These values also would in- 
dicate that the cooling break, v c , was located at frequencies 
lower than the R-band ~ 8 hours after the gamma-ray event. 
This would make from GRB 00101 1 a very interesting system, 
since in most cases v c is higher than the optical frequencies at 
early times. 

Images taken 7 months after the burst reveal an elongated 
object with R= 25.38 ± 0.25 fully consistent with the OA po- 
sition, likely the host galaxy of GRB 00101 1. 
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Fig. 2. Upper left panel: The panel shows an R— Ks versus J— Ks colour-colour diagram for the objects inside the GRB 00101 1 
error box (open diamonds). The star represents the colours of the GRB 001011 optical/near-IR counterpart. The errors of the 
counterpart colours are smaller than the size of the star. The shaded region shows the OAs colour-colour space locus determined 
by our selection criteria R— Ks > 2.0, J— Ks > 1.2. As can be seen, the colours of the GRB 001011 counterpart are consistent 
with the shaded area and inconsistent with almost all the other objects within the error circle. Only two non-transient objects in 
the error box have colours consistent with the shaded area. They are very likely quasars or compact galaxies (see Sect. 3.4 ). Upper 
right panel: The figure shows an R— Ks vs J— Ks colour-colour diagram for the 10097 2MASS+USNO sources (dots) found in a 
±5° window around GRB 00101 1. As it is shown, the counterpart (star) exhibits colours inconsistent at least at a 3a level with 
the objects in the 2MASS+USNO catalogue (see Sect. 



3.1 



for information on the iso-density contour levels). As can be seen, the 
2MASS+USNO catalogue traces most of the colours shown in the upper left panel. The seven open squares labelled from 1 to 7 
represent the colours of afterglows with measured R— Ks and J— Ks colours to date. The additional square labelled with number 
8 represents GRB 000131 which was not detected in the J-band (Andersen et al. 2000). As can be seen, at least seven of the eight 
open squares and the star are consistent with the shaded area. The eighth open square (GRB 000131) is likely consistent with this 
region too. The solid straight line inside the shaded region represents the colour trace of a pure power-law SED with (3 ranging 
from 0.6 to 1.5. When Lyman-a blanketing is taken into account the pure power-law SED (straight solid line) is shifted upwards 
(three dashed lines). From the bottom to the top, the three dashed lines show the colours of a power-law SED when the Lyman-a 
blanketing is considered for an afterglow at z = 5.0, 5.5 and 6.0, respectively. As can be seen, the higher the redshift the more 
distant the power-law SED from the stellar trace is. Lower left panel: The figure displays synthetic colour-colour traces for stars 
(filled circles), galaxies (dashed curves) and quasars (dotted curves). The shaded polygonal background region represents the loci 



occupied by the asteroid samples of McCheyne, Eaton & Meadows (1985 1 and Gaffey et al. (1993). As can be seen, the bunch of 



